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Abstract 

A supersymmetric standard model with heavier scalar supersymmetric particles 
has many attractive features. If the scalar mass scale is O(10 — 10 4 ) TeV, the 
standard model like Higgs boson with mass around 125 GeV, which is strongly 
favored by the LHC experiment, can be realized. However, in this scenario the 
scalar particles are too heavy to be produced at the LHC. In addition, if the scalar 
mass is much less than 0(1O 4 ) TeV, the lifetime of the gluino is too short to be 
detected. Therefore, it is hard to probe the scalar particles at a collider. However, 
a detailed study of the gluino decay reveals that two body decay of the gluino carries 
important information on the scalar scale. In this paper, we propose a test of this 
scenario by measuring the decay pattern of the gluino at the LHC. 



1 Introduction 



A supersymmetric (SUSY) standard model (SSM) is the most promising model beyond the 
standard model (SM). At the Large Hadron Collider (LHC) experiment, the Higgs boson 
search of the ATLAS and CMS collaborations reveals the existence of the Higgs boson- 
like particle with mass around 125 GeV pQ El |3]. However, the minimal SSM (MSSM) 
naturally predicts relatively light mass of the Higgs boson, and thus the parameter region 
where the Higgs mass is heavy is strongly constrained. To realize the Higgs mass about 
125 GeV, a heavy scalar top quark (stop) and/or large A-term of the stop are required, 
in the case of the MSSM 

Among them, heavy scalar scenarios are interesting from many viewpoints. The Higgs 
mass around 125 GeV indicates the supersymmetric scalar mass scale of O(10 — 10 4 ) TeV 
[6j [7J EJ |9]. Such a heavy scalar particles can greatly relax the notorious constraints of 
the SUSY flavor and CP problem and the cosmological gravitino problem as discussed 
in the context of the split SUSY pm HU [12], the anomaly mediated SUSY breaking 
(AMSB) models [T3| [T4] and the spread SUSY [15]. From a theoretical point of view, 
it is more likely that the scalars have larger mass compared to the fermions, such as 
gauginos and Higgsinos. This is because the mass of fermions can be protected from a 
large radiative corrections thanks to some symmetry and, on the other hand, that of the 
scalar particle is not. Therefore the MSSM with heavy scalar particles is favored from 
both phenomenological and theoretical viewpoints. 

In order to test this scenario, an estimation of the scalar mass scale, especially for 
the stop mass, is quite important, since it strongly connects to the Higgs mass. This 
also provides insights into the low-energy precision measurements, such as flavor and 
CP physics. However, the scalar mass is too heavy to be produced at the colliders in 
this scenario, even if the gauginos are light enough to be produced. Therefore we need 
indirect information to study the scalar sector. For example such heavy particles can 
be relevant in the early universe, and, in some case, this mass scale can be determined 
by a gravitational wave [16]. A precise measurement of gaugino-Higgsino mixings at the 
International Linear Collider (ILC) may enable us to estimate the scalar mass scale by 
using renormalization group (RG) equations [17] ■ At the LHC, if the scalar mass scale 
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is much larger than about 10 3-4 TeV, an R-hadron or a displaced vertex of a long-lived 
gluino is expected [TUl EL71 HE] • The lifetime of the gluino is very sensitive to the scalar 
mass scale, and hence measurement of the lifetime of the gluino will be a probe of such a 
heavy scalar mass. However, this method cannot be applicable in the case that the scalar 
mass scale is much smaller than O(10 4 ) TeV, since the decay length is too short to be 
detected at a collider experiment. 

In this paper, we study a possibility to test the scalar mass scale in the heavy scalar 
scenarios through the gluino decay patterns. In the case of the short-lived gluino, what 
we can measure is basically branching fractions of the gluino. If we study only the gluino 
three body decayJ3 it is difficult to extract the information of the scalar mass scale. On the 
other hand, the two body decay of the gluino is sensitive to the scalar mass scale, as we 
will see in the next section. Actually, the Higgs mass around 125 GeV implies a possibility 
of a large size of two body decay of the gluino. We discuss that the branching fraction of 
9 ~~ ^ 9X° carries important information of the scalar mass scale, and we propose a simple 
method to measure it. This method relies on a naive counting of the events which are 
classified with their jet, lepton and 6-jet multiplicities and kinematical behaviors of jets 
and missing energy. In order to define signal regions, we only use commonly-used cuts 
adopted in the study of the ATLAS and CMS collaborations. We show that the branching 
fraction can be determined to accuracy of 10 % if the (9(10 4 ) gluinos are produced at the 
LHC. 

This paper is organized as follows. We first review the gluino decay in the case of the 
heavy scalars in Section 2. We propose a basic strategy to discriminate the two body 
decay of the gluino and demonstrate it for several models in Section 3. The conclusion 
and discussion are given in Section 4. 

2 Heavy scalar scenario and gluino decay 

In this section, we briefly review the gluino decays in the heavy scalar scenario. In this 
scenario, we can neglect the left-right mixing of the squarks, because the mixing angles 

1 Previously, the three body decay of the gluino in the heavy scalar scenario is studied in the different 
context. See, for examples, Refs.jl9. . 
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(a) Three body decay (b) Two body decay 

Figure 1: Examples of diagrams relevant for the gluino decay. 



are suppressed with quark masses and are negligible in the present situation. 

2.1 Leading order calculation 

Three body decay (g — > qqx) 

Squark exchanging diagrams at tree level induce three body decays of the gluino (Fig. 
[H-(a)). We assume the squark sector respects flavor symmetries for simplicity. For 
massless quark limit, we can get an analytical expression of the gluino partial decay 
width, e.g., 

A<) + ^ 9 (<)}, a, 

where f(x) = 1 — 8s — 12x 2 logx + 8a; — x 4 and g(x) = l + 9x + 6a; logs — 9x 2 + 6x 2 logx — x 3 . 
As we can see from the above expression, the decay widths simply damp as m^ 4 in the 
heavy squark limit. 

Two body decay (g — > gx°) 

The gluino can decay into a neutralino and a gluon via quark-squark loop diagrams 
(Fig. [T]-(b)). First, one can easily find a branching fraction to a Wino (W) is strongly 
suppressed, because such a decay mode is induced by dimension 7 effective interactions, 

2 In the heavy scalar scenario, the flavor violation in the squark sector is less constrained, because the 
dangerous flavor-changing-neutral-current (FCNC) processes are suppressed by the large squark masses. 
However, in spite of the heavy squark, the gluino three body decay can be a good probe of the flavor 
violation in the squark sector, because three body decay is a direct result of the four-Fermi interaction. 
We will study this possibility in more detail elsewhere [20]. 
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such as, 

777 — 

£ cff . ~ (H^r a H)W a a^~gG, u . (2) 

Then, g — > gW is mainly induced by the mixing with other electroweak(EW)-inos. 

As for the gluino decay into a Bino (B) or a Higgsino (h), the decay rates T(g — > gB) 
and T(g — > gh) are given by, 





40 967T m~ \ mf J m \ R 

Here, y t is Yukawa coupling constant in the MSSM. The index q in Eq. fl3]) runs through 
u, d, s, c, b and t, and Y q is the hypercharge of the quarks. In Eq. (j4]), we neglect contri- 
butions other than top-stop loops, because it has a dominant contribution. Eq. fll]) shows 
r(<7 — > gh) has an enhancement factor mj/m?(log(m?/m^)) 2 , although T(g — > g5) and 
r(p ggx) simply damp as w^ 4 . Then, a rough behavior of the ratio of T(g — > gh) to 

^4«i(, og <) 2 . (5) 

r(0 -> m | V m tJ 

Therefore, we can expect the branching fraction of the gluino two body decay is sig- 
nificantly enhanced in the heavy scalar scenario. The logarithmic term in Eq. (pE} can 
be understood as an effect of the RG from the view point of the effective theory after 
integration of the heavy squarks [21] . 

2.2 RG-improved calculation 

In the present model, the squark mass scale is much larger than the gluino mass and thus 
the above leading order calculation is not suitable, especially for the gluino two body decay. 
To estimate the gluino decay width more precisely, we integrate out the squarks and use 
the effective Lagrangian which describes the dynamics of the SM particles together with 
Higgsinos and gauginos. In this setup, the gluino decay is caused by non-renormalizable 
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interactions which are induced by integrating out the squarks. We discuss the Wilson 
coefficients of dimension 5 and 6 operators which govern the gluino two body decay at 1 
loop level. Such operators are given by, 

C e s. = C?Q* + (ci0§ + C*Q* + h.c.) . (6) 

Here, Q's are operators defined below: 

Qf = Ba^g G^, (7) 

Q§ = [Ha^P R g a ) {q^ ] a, u T a P R t) , (8) 

Ql = Ha^P R ~ghG, u . (9) 



We obey the notation of Ref. [21]. We have to match the low energy effective theory to 
MSSM at some matching scale m, which is set to a typical scale of the squark masses. 
The Wilson coefficients at the matching scale rh are given by, 

- it^K-^>(£^-^f)> ("» 
eft*) = jjf—A -4- + 4- ], (id 

C * {Th) = vtl' R + (12) 

32\/27r 2 smp \ m_ C3) m- 

\ Or l R 



Here, q in Eq. (fTUjl runs through u, d, s, c, b and i. y t is a coupling constant for top-top- 
Higgs coupling like the SM. We can see that Cf(fh) and C^{m) are suppressed by a loop 
factor compared to C|^(m). Below the scale fh, the running of the Wilson coefficients is 
described by the following equations: 

16vrV^-Cf = -Ug 2 s C*, (13) 




fg 2 + \y 2 t 2g s y t 
Ag s y t -Ug 2 s + 3y 2 




(14) 



Eq. (|14p shows, in spite of the suppression by a loop factor, C§ is significantly enhanced 
at low scale because of mixing between Q2 and Qf. On the other hand, Cj is not 
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enhanced because does not mix with four Fermi operators at 1 loop level. § A 
neutralino mass eigenstate is a mixing of the Bino, Wino and Higgsino. Then, the decay 
width of the gluino two body decay can be written by the mixing matrix iVj.,- and the 
Wilson coefficients at low scale. □ The decay width of the gluino two body decay is given 



by, 



(rria — m^o) / 3_vv° n \ - 
-> tfy) = o_ 3 ' > (15) 

2 



27rm~ 



cJ^'C/i) = CfiriNa + cii^Nuv + ^ciiriNuvlog 1 ^-. (16) 

Here, v ~ 174 GeV is the SM Higgs VEV. We have checked the above expressions are 
consistent with Eq. ([3]), (j4j) at the leading order. 

2.3 Numerical evaluation of the branching fractions 

In closing of this section, we show some numerical evaluations of the branching fractions. 
We use the RG- improved calculation by using the RG equations given in Refs. [7J [21]. 
Fig. [2] shows Br(<7 — » decay length ct~ 9 and the lightest Higgs mass rrth in the 

heavy scalar scenario. We can see an enhancement of the branching fraction of the two 
body decay as the stop mass increases, and there exist parameter regions consistent with 
rrih — 125 GeV. An important observation is that we can expect more than O(10) 
% branching fraction of the two body decay when tan/3 = 0(1). The partial decay 
width of the two body decays is almost determined by the stop mass scale, although the 
total decay width is sensitive to the first and second generation squark mass and so on. 
Hence, to extract the stop mass scale, Tl g /t,b = r(p — > gx)/F(g —> tt/tb/bb + x) is more 



3 This can be seen by explicit diagram calculation. By using Fierz identity, the four Fermi operators 
which contain a gluino and a Bino can be written by, 

qLh = (Bri 5 m^T a p L . Rq ). 

At 1 loop level, only a quark loop diagram is a candidate to give mixing of Q^ L and Qf. However, the 
quark loop contribution is proportional to the correlation function of a quark current, then, its transverse 
part is proportional to g 2 , where q is momentum of a gluon. Therefore, when an external gluon is on-shell, 
its contribution is vanished. This behavior is not guaranteed by any symmetry, therefore, Q^ L and Cf 
have mixing at 2 loop level. 

4 In this paper, we neglect the radiative corrections on the neutralino and chargino mixings from the 
scalar particles for simplicity. 
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suitable parameter than Br(g — > gx)- In Fig. [31 we show a scatter plot of the stop mass 
m sto P = ^ m ql m t R a t the stop mass scale versus Tl g /t,b- In this figure, we take m g = 1 TeV 
and other parameters are chosen from a parameter region: 

• 100 GeV < Afi, M 2 , /i < 400 GeV. 

• 2 < tan/3 < 10. 

• 3 TeV < m stop < 10 5 TeV. 

• 0.5m^ R < 77^1,2,77^3,771^ < 2m^ at the Grand Unified Theory (GUT) scale. 

Fig. [3] shows that Tl g /t,b which is mainly determined by m stop , and the information about 
the EW-ino masses enable us to determine m stop more precisely. For larger fi, the two 
body decay is relatively enhanced because the sizable branching fraction of three body 
decay into Higgsino and quarks becomes smaller. Fig. 0] shows how Ti g /t,b depends on 
the gluino mass rrig, and this figure shows the light gluino is advantageous to enhance 
H g /t,b- As we can see from Eq. (j3J or Eqs. (1TS| [TBI) , the enhanced part in T(g ^ gx) 
needs chirality flip by the top quark mass. Therefore, T(g — > gx) has a suppression factor 
m t/ m $- Finally, we show a sample point of the gluino decay table and mass spectrum in 
Table. OQ 



Table 1: We show the sample point of the gluino decay table (left) and mass spectrum 
(right). We take M 1 = 150 GeV, M 2 = 300 GeV, /i = 230 GeV, tan/3 = 2.5 and 

0.123, 



niq = 2 x 10° GeV. We assume universal squark pole masses. Br(g — > gx 



n g/ttb = 0.459. 
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Figure 2: Br(2 body) (solid red lines) and ct~ 9 (chain blue lines) as a function of tan/3 
and m st0 p. In dark gray shaded region, 124 GeV < rrih < 126 GeV. The light gray region 
shows the uncertainty of rrih due to the present error on the top quark mass. At the stop 
mass scale, we take first and second generation squark mass as m stop (left) and 2m stop 
(right). We take M x = 150 GeV, M 2 = 300 GeV, M 3 = 1000 GeV, li = 230 GeV and 

and a s {m z ) = 0.1184 [23]. 



173.2 ±0.9 GeV (22 




Figure 3: Third generation squark mass versus TZg/t,b = T(g — >■ gx)/V{g — > tt/tb/bb + x)- 
On the black line, we take M x = 150 GeV, M 2 = 300 GeV, /i = 230 GeV and tan/3 = 2, 
and the gray lines show the variation if we change the parameters by 20 %. We take 
uniform distributions for Mi,M2,tan/3, m^ li2 /mj , m%/mj and m~ /m? at the GUT 
scale, while we take In m st0 p to be uniformly distributed for m st0 p = y/ m ql m t R - 
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Figure 4: Tl g /t,b = T(<? — > gx)/Y(g — > tt/tb/bb + x) as a function of the gluino mass. We 
take the stop mass as 10 4 , 10 6 and 10 8 GeV. For other parameters, we take M\ = 150 GeV, 
M 2 = 300 GeV, \i = 230 GeV and tan/3 = 2. We take universal squark masses at the 
squark mass scale. 



3 Study of the gluino decay at the LHC 

As we discussed above, the branching fraction of the gluino two body decay carries im- 
portant information on the scalar mass scale. In addition, the recent Higgs observation 
suggests sizable branching fraction of the gluino two body decay. In this section, we dis- 
cuss the measurement of the gluino branching fractions at the LHC. For this purpose, 
we adopt a simple method, using the difference of the event topology depending on the 
gluino branching fraction. Then we demonstrate this method for some models. 

3.1 Basic Strategy 

First we compare the three body decay (g — > q^x) mediated with third-family squarks 
to the two body decay into a gluon (g — > gx)- For the former case, the gluino decays 
into third-family quarks, b and t, if small flavor violations assumed. These 6-jets and 
leptons originated from the top quark decay will leave distinct signatures from the latter 
two body decay case. Therefore, it is relatively easy to discriminate the two body decay 
from the three body decay into the third-family quarks. 
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Then, let us discuss the discrimination between the two body decay and three body 
decay into the first and second-family quarks. Unlike the third family case, jet flavors and 
leptons are not necessarily help for the discrimination. However, the behavior of the jets 
Pr and missing energy will be distinguished because of the difference of kinematics of the 
two body and three body decay. Naively, the three body decay case is expected to have 
more number of jets and each jet has less energetic compared to the two body decay case. 
To extract the difference, we define event cuts according to jet multiplicity and study 
the jet behavior together with E T:miss in detail. These cuts enhance the difference of two 
decay patterns and also take a role of reduction of the SM background as well. 

Another barometer will come from difference of the neutralino in the gluino decay 
product. As for the decay g — > gx°, as discussed before, the final state x° have large 
H u component, since the two body decays requires the large Yukawa coupling of the top 
quark and the neutralino. On the other hand, the three body decay mediated with the 
first and second family squarks G?m,2 and Um t 2 responds to the gaugino components 
B and W . Therefore the species of final state neutralino from the two body and three 
body decay can be used as an information to distinguish the decay patterns of the gluino. 

Using these differences, we can estimate the gluino branching fraction. We adopt the 
following method. One of the basic strategy is owe to a simple number counting of the 
leptons and 6-jets. 

We separate signal regions corresponding to the number of the leptons Ne, r-jets N T 
and 6-jets N b in each event. Here we require that the P T 's of the lepton, r and b are larger 
than 30 GeV, 50 GeV and 50 GeV, respectively. In addition, as for the lepton numbers, 
we consider set of N% = {0, 1, SS2£, OS2£, > 3}, where SS2£ represents same sign 2 leptons 
and OS2£ opposite sign. As for the r and 6-jets, N T = {0, > 1}, iVj, = {0, 1, 2, > 3}. 

We also count the number of jets Nj with Pt > 100 GeV. To characterize kinematical 
nature of the event, we further consider two conditions referred to as A^-Jl and A^-J2 
for each jet multiplicity and each event is tagged whether it satisfies these conditions or 
not. After all, each event can be labeled with a set of (N b , N t) N T , Nj, A^-Jl, Nj-32). We 
always require Pr,miss > 300 GeV as a basic cut. 

For multi-6-jets and leptons modes, this basic cut and requirement of multi 100 GeV 
Pr jets reduce the SM background effectively. However, in the case of less 6-jets and 
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leptons events, the basic cut is not adequate for the SM background reduction and also 
is difficult to discriminate two body decay and three body decay. Then the conditions 
Nj-Jl and Nj-J2 play an important role. They are defined as follows. 

2 jets mode Nj = 2 

A0(ji i2 ,-&r,miss) > 0.4 rad and 

2J1: P jliT > 500 GeV, P j2 , T > 250 GeV, E T , miss > 600 GeV, M T2 > 650 GeV 
2J2: P j2 , T > 150 GeV, M eff>2 > 1200 GeV, E T<miss > max(0.25M effi2 , 400 GeV). 

3 jets mode Nj = 3 

A0(ji i2;3 ,E Timiss ) > 0.4 rad and 

3J1: P jl;T > 350 GeV, P j2 , T > 300 GeV, M eff>3 > 1200 GeV 
3J2: P j2T > 150 GeV, M cffi3 > 1500 GeV, E T<niiss > 0.35M effi3 . 

4 jets mode Nj = 4 

A0(ji, 2 , 3 ,4, ^T,miss) > 0.4 rad and 

4J1: P AT > 250 GeV, P i4 ,T > 200 GeV, M e $ A > 1400 GeV and E TMss > 0.1M effi4 . 
4J2: P j2T > 150 GeV, M cfM > 1400 GeV, E TMss > max(0.35M cffi4 , 600 GeV). 

5 jets mode Nj = 5 

A0(ji, 2 , 3 ,4, E TMss ) > 0.4 rad and 

5J1: P jltT > 400 GeV, P AT > 150 GeV, M eff>5 > 1400 GeV, 

P T ,miss > max(0.25M effi5 ,700 GeV). 
5J2: Pj 2T > 300 GeV, P j3>T > 150 GeV, P T ,miss > 0.1M effi5 . 

6 and more jets mode Nj > 6 

A0(ji, 2 ,3,4, pT,miss) > 0.4 rad and 

6J1: P jhT > 200 GeV, P i6 , T > 150 GeV, M effiall > 1400 GeV. 
6J2: P j3T > 150 GeV, M eSM i > 1400 GeV. 

Here P jrijT is the n-th highest P T of jets and M eSti = Y,i< n <i P jn,T + E TMss + Si ep tons ^,t- 
Note that the cuts imposed above are pretty normal, and actually similar cuts are already 
used in the current analysis [25]. By using the number of events after each cut, we can 
confirm the above qualitative arguments in a quantitative form. 



12 



Then let us consider the determination of the branching fraction. We use the maximum 
likelihood estimation (MLE). We consider the following likelihood C, 

C= J] Prob(iV mode |{ J B l ,^}), (17) 

modes 

where Bi is a set of branching fractions of the gluino, a- gg is the gluino production cross 
section and N mo ^ e is the number of the events characterized with (iV 6 , Np, N T , Nj, Nj-Jl, 
Nj-J2). We assume Prob is a convolution of Poisson distribution from the event number 
and Gaussian distribution from systematic uncertainties of the acceptance of the SUSY 
signals and the number of the SM backgrounds. In this study, we fix the parameter such 
as the gluino mass except for the gluino branching fractions and the cross section. As 
for the systematic uncertainties of the number of the SM 5b, we assume 5b = 20 %. For 
acceptances of the SUSY events (5s), we set 5s = 20 %. In the above set-up, we consider 
the minimization of — 21n(£) by varying {Bi,a 9 g} and estimate best fitting parameters 
and their errors as the usual MLE method. 

In the following, we show the result of Monte-Carlo simulation. To estimate the SM 
background, we have used the programs MC@NLO [26J (for tt, WW, WZ and ZZ), Alpgen 
[27] (for Wj, Zj and W/Z + bb/tt) and Pythia 6 [28] (for QCD jets). As for the detector 
simulation, we used the AcerDet [21]. For the isolation criteria for leptons and photons, 
we used the default setting of the AcerDet. In our study, we are interested in the signals 
with leptons, r jets and 6-jets. The detection efficiency and misidentification of these 
particles are important in estimating both the signal events and background events. In 
the following simulation, we include the misidentification and fake rates of the leptons, r, 
fo-jets and light jets, following Ref. [21]. As for the MSSM mass spectrum, we have used 
the program ISAJET 7.82 |30j modified for the scalar decoupling scenario. The SUSY 
events are also generated with the program Pythia [28J . The NLO production cross section 
of the gluino production is estimated with the program Prospino2 [3T] . 

3.2 Examples 

We apply the above method on three models as examples, and see how accurate we can 
determine the branching fraction of the gluino two body decay which is sensitive to the 
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mass scale of squarks. We consider a simplified model which contains only a gluino, a 
neutralino and a chargino, the split SUSY model and the anomaly mediation model. 

3.2.1 Simplified Model 
Setup 

Let us consider only three particles are relevant: gluino g (fOOO GeV), a light neutralino 
(100 GeV) x° an d a chargino X s1 slightly heavier than the neutralino (by 100 MeV). This 
simplified set-up is very useful to see the how well the above method works. We set the 
chargino decay x — > eux°, whose decay products are irrelevant for the collider study. 
The collider signature can be characterized by the branching fractions of the gluino: 
g — y gx° , ttx° , tbx , bbx° an d qqx°- Thus, we vary branching fractions of these decay 
modes, Bi, as well as the production cross section of the gluino-pair, o~ g ~ g - 

Measurement of the Branching Fractions 

We apply the strategy discussed in Section 3.1 to this simplified model, and determine a 
set of {Bi, (jgg} based on the MLE. To see the difference of the event topology intuitively, 
it is useful to use the significance variable Z, which represents the deviation from the 
SM background and roughly given by ratio of the number of SUSY events to that of the 
fluctuations of the SM background [32] . 

In Fig. [51 we show the significance variables Z's of some selected cuts which is im- 
portant for the parameter estimation. Here we assume an integrated luminosity of 5 fb _1 
at -y/i = 14 TeV LHC run. As expected, the gluino decay into the third-family quarks 
provides very distinct signatures from the gluino two body decay. As for the difference 
of decay into qq and g, we can see the clear difference when we compare signal regions of 
higher and lower jet multiplicities. The the decay mode of g — > qqx° shows large Z for 
large jet multiplicity modes. 

In Fig. |6l we show 1 and 2-a contour lines of expected parameter estimation for 
combination of B{g — > gx) and B(g — > qqx) with an integrated luminosity 10 fb _1 at 
y/s = 14 TeV. Here we assume that B(g — > gx) — 0.3, B(g — > tlx) = 0.7 and others 
as fiducial values. 
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Figure 5: Significance variable Z for each mode in the simplified model. In the row of 
Jets, iVJ represents the number of the jets with Pt > 100 GeV and iVJl,2 represent the 
additional conditions in Sec 3.1. 
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Figure 6: 1-a and 2-a contour lines of the estimations of B(g — > gx) and B(g — > ttx) at 
the sample point in the simplified model. The point x represents the fiducial value. 



There is a direction, B(g — > gx) + B(g — > qqx) — 0.3, in which it is slightly difficult 
to determine each branching fraction since both of them do not generate 6-jets or leptons. 

In Table [2} we show the expected accuracy of determination of the branching fraction 
of two body decay for some sets of the parameters. In the simplified model, we can see 
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the branching fraction of two body decay can be determined to accuracy of 0.02 - 0.05. 
Table 2: Expected accuracy of two body branching fraction for the simplified model. 



Fiducial Values 


£ [fb- 1 ] 


lo" estimation 


(B(g),B(tt)) — (0.3,0.7) 


10 


(0.26, 0.32) 


(B(g),B(tt)) — (0.3,0.7) 


25 


(0.27, 0.31) 


(B(g),B(qq)) = (0.2,0.8) 


10 


(0.15, 0.24) 


(B(g),B(qq)) = (0.2,0.8) 


25 


(0.17, 0.23) 



3.2.2 Split SUSY 
Setup 

Then let us consider a more realistic example. Here, we discuss the scalar decoupling 
scenario. At the low energy, there are the gauginos (B, W and g) and the Higgsinos 
(HujHd), which are relevant for the LHC study. We assume the Higgsino mass is the 
same order as the Bino and Wino masses, and then, EW-inos and Higgsinos are mixed 
with each other and form the mass-eigenstates: neutralinos and charginosl^l We assume 
that the gauginos masses and the mixing parameters are well-known. This assumption 
would be justified by considering that the neutralino can be the dark matter. For example, 
the abundance, the cross section to nucleon and mass and/or cosmic ray signatures at 
future experiments provide information on the EW-ino sector. Kinematic information 
at the LHC also gives a large number of clues to determine this SUSY EW-ino sector 
[34"1 [35] . Or, in an extreme argument, an e + e~ linear collider such as the ILC can provide 
the precise information on the neutralino and chargino sector [17]. As for the neutralinos 
and charginos decay, the contribution from the heavy scalars are negligible when the scalar 
mass are much heavier than TeV scale. In this case, the branching fraction of the gluino 
two body decay is a few ten percent, if universal scalar masses assumed. 

5 To be precise, the mixing parameters in the split model arc different from ones of the usual MSSM 
scenario, because of the radiative corrections due to the large mass hierarchy between the squarks and 
the gauginos. However, we neglect this difference for simplicity. 

6 In this set-up, the dark matter abundance fth 2 is around 0.1. If the scalar scale is around 1000 TeV, 
the SM-like Higgs mass is around 125 GeV. 
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Then we consider the gluino decay. As discussed in the previous section, the gluino 
decays are induced by the dimension 5 and 6 operators. However, there are a lot of the 
operators which induce the gluino decay, and then, it makes the analysis very complicated. 
To make the analysis easier, we take two simplifications. First simplification is a limitation 
of the operators. We only introduce the following seven operators: 

= c^B^mi'PRq) + c tRtR (BYi 5 g)(iYPRt) + c bRbR (BYi 5 g)(h"P R b) 
+ c qLqL (WYi 5 ~g)(qYPLq) + c Q 3 Q3 {Wr 1 5 g)(Q 3 rp L Q 3 ) 

+ C Q3tR (H u g)(tP L Q 3 ) + C g G, u H° u a^g, (18) 

where q collectively represents first and second generation quarks. Given the seven C"s 
in Eq. (JTBl) . Mi, M 2 , /x, tan/3 and m§, we can calculate the decay width of each decay 
mode. For example, let us consider T(g — > tix )- This decay width is determined by 

c t R t R , Cqsqs and C Q 3 tR , such as, 

T(g ->• tix ) = C 2 tRtR A tRtR + C 2 Q3Q3 A Q 3 Q3 + C 2 Q3tR A Q3tR (19) 

+ Ct R t R CQ3Q'iB tRtRj Q3Q3 + C tRtR C Q 3 tR B tRtR) Q3 tR + CQ3Q3CQ3 tR BQ3Q3 iQ 3 tR , 

where .A's and £>'s are functions of the masses of the quarks, neutralinos and charginos. 
The first line of the above equation shows a linear combination of the contribution in case 
that only one operator is switched on and the others off. The second line shows interference 
effects between the operators due to the top quark mass and the EW-ino mixings. The 
first line can be calculated easily by a linear combination of the contribution from each 
operators, although the second line can not. Hence, we take a second simplification. We 
drop the interference terms between different operators, i.e., the second line in Eq ffT9]) . 
Such a simplification is justified if the gluino mass is much larger than the top quark 
mass and the EW-ino mixing terms, and it makes the gluino width into a simple linear 
combination of the decay width in case that only one operator is introduced. Similarly, 
the branching fraction can be written by a linear combination. For example, if we denote 
the branching fractions Brj(^ — > ttx°) in case that only Cj is switched on, the general 
branching fractions is given by, 

Br(^ tix ) = BiBl ^9 ->• ttx°), (20) 
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corresponds to the case that the coefficient Ci in Eq. f|T8l) is dominated. 



where {B,{\ = { B% B^, fl^, B' QaQa , 2%^, S^, Bj^} with = 1. B.'s have 

common value for all decay modes. B^s are weights in a linear combination, and they are 
determined as the functions of the Wilson coefficients Cj . B 2 g corresponds to the branching 
fraction of the two body decay of the gluino, because the only operator C g G^H^a^g 
in Eq. (ITS]) is able to induce the two body decay. From now on, we tackle with an 
optimization problem of Bi instead of C». Here, we set nig = 1000 GeV, M<i = 300 GeV, 
Mi = 150 GeV, n = 230 GeV and tan/3 = 2 as a sample point. 

Measurement of the Branching Fractions 

Then let us move to the parameter estimation. Under the above assumptions, we esti- 
mated accuracy in determining the gluino branching fractions at the LHC, using the MLE. 
In Fig. [7J we show significance variables Z's for some cuts for C = 5 fb _1 at y/s = 14 
TeV. As in the case of the simplified model, the two body decay mode of the gluino has a 
smaller Z for high jet multiplicity modes. As expected, the three body decays into third 
family quarks are distinguished by the number of 6-jets and/or leptons. 

In Fig. |H1 we show 1 and 2-a contour lines of expected parameter estimation for B 2 g 
and + B^ . Here we assume that B(g — > gx) = 0.3, B(g — > ttx) = 0.7 and 
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X 

0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Figure 8: 1-cr and 2-cr contour lines of the estimations at the sample point in the split 
SUSY model in which B(g — > gx) = 0.3 and B(g — > tlx) — 0-7. The point x represents 
the fiducial value. 

others and that an integrated luminosity 25 fb _1 at y'i = 14 TeV. As in the case of the 
simplified model, we can see there is a direction, B 2 g + B% + B^ lQl ~ 0.3, in which it is 
slightly difficult to determine each branching fraction. 

In Tab. El we show the expected accuracy of determination of the branching fraction 
of two body decay. We also show the case that some of the explanatory variables are 
constrained. As expected, the gluino three body decay into the first and second-family 
quarks mimics the signals of the gluino two body decay. With ad hoc assumption of the 
absence of such three body decay, a large improvement of the parameter estimation is 
possible. Once we get the information about the gluino decay branching fraction, we can 
get the stop mass from 1Zg/t,b defined in the previous section. For example, if we take 
Mi = 150 GeV, M 2 = 300 GeV, /i = 230 GeV and tan/3 = 2, we can get the stop mass 
from the black line in Fig. [3j From Tab. |3l we can get an uncertainty STZ g / t ,b is roughly 
given by 0.05 at 50 fb~ . This uncertainty corresponds to the stop mass uncertainty 
51ogm stop ~ 0.7, then, roughly, we can get the stop mass with the uncertainty of the 
factor 2. 

Let us comment on the uncertainties from the EW-ino sector. In this study, we assume 
that the parameters of this sector are precisely determined. However, this is not a crucial 
assumption at the determination of the branching fractions. The present method mainly 
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Table 3: Expected accuracy of determination of the two body branching fraction of the 
gluino in the split model. 



Fiducial 


£ [fb- 1 ] 


1-a estimation 


Condition 


(Bt,B* KQn ) = (0.3, 0.7) 


25 


(0.23, 0.32) 




[Bl,Bl Q J = (0.2, 0.8) 


25 


(0.14, 0.22) 




{BlBl Q? ) = (0.3,0.7) 


50 


(0.25, 0.31) 




(5^?^) = (0.2, 0.8) 


50 


(0.16, 0.21) 




(^,5^) = (0.3, 0.7) 


25 


(0.28, 0.32) 


B 3 = B 3 =0 

QBQR QLQL 


(£^ Q ,^,aJ = (0.2,0.75,0.05) 


25 


(0.14, 0.24) 




(£^Q.^LJ = (0.2,0.75,0.05) 


50 


(0.16, 0.23) 





relies on the the third family quarks and high Pt jets from the gluino decay and not on 
the detail of the neutralino and chargino sector so much. Actually, we have checked that 
0(10)% uncertainties of the parameter of the EW-ino sectors do not change the result 
significantly. On the other hand, at the determination of the stop mass from TZ g u^, Fig. |3] 
shows us that (9(10)% uncertainty about the parameters in the EW-ino sector introduces 
the factor 2-3 uncertainty about the stop mass. Then, we do not need the detail of the 
EW-ino parameters to get the order of the stop mass, however, the precise determination 
of the stop mass requires both the EW-ino sector information and more large luminosity. 
The present method is robust but crude, and there can be room for improvement. When 
addressing this improvement, we must deal with uncertainties from EW-ino sectors, which 
is out of the scope of this paper. 

3.2.3 Anomaly Mediation 
Setup 

Another important scenario of the heavy scalar scenario is anomaly mediation models. 
In this model, the gaugino masses are generated from heavy gravitino mass thorough a 
conformal anomaly. The gaugino mass from the gravity mediation can be suppressed, if 
appropriate symmetry assumed. In this case, the Wino is the lightest supersymmetric 
particle (LSP). On the other hand, it is difficult to forbid the scalar mass generation from 
the gravity mediation effects by using some symmetry. Therefore a natural prediction of 
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the AMSB model is heavier scalar mass around the gravitino mass (0(100) x m gaug i no ) or 
more. Unlike the split SUSY case, we assume that the Higgsino mass is as heavy as the 
scalar particles. 

In the case of the AMSB model, the charged Wino can be long-lived because of small 
mass difference between charged and neutralino Winos. This charged track can be de- 
tected at the LHC if a large number of SUSY particles have adequately large cross section, 
which is provides an important test of the AMSB model [361 [371 [38] . However, the decay 
length of the charged Wino is typically small (9(1 — 10) cm and dominant signals of the 
AMSB models are just typical SUSY signals, i.e, large Pt jets and large missing energy. 
Here we focus on this dominant AMSB signals. 

Let us consider the gluino decay in the AMSB model. The gluino two body decay 
into a neutralino and a gluon mainly comes from the Higgsino component. Thus in this 
AMSB case, decay g — > g\o is significantly suppressed. Therefore the AMSB predicts 
absence of the gluino two body decay and its confirmation provides an important test 
for the AMSB model. Here we discuss how well the gluino two body decay can be 
constrained. The gluino decay can be essentially parametrised by following branching 
fractions: g -> gW ,gI3 ,tiW ,bbW ,tbW ± ,qqW,tiB ,bbB°,qqB . As a sample point, 
we assume rrig = 1000 GeV, m^ = 450 GeV and = 150 GeV. 

Measurement of the Branching Fractions 

Using the method described previously, we constrain the gluino two body branching frac- 
tion. In Fig. [HI we show significance variables Z's for some cuts for C = 5 fb _1 at y/s = 14 
TeV. As in the case of the simplified model, gluino two body mode have smaller Z for 
large jets multiplicities modes. In Table HI we show the expected upper bound on the 
two body branching fraction. In this case, we can see that the gluino decay width well 
constrained, which gives good test for the AMSB scenario. 

4 Conclusion and Discussion 

The gluino two body decay into a gluon is a quite important feature of the scalar decou- 
pling scenario. As the scalar mass scale gets larger, the gluino two body decay fractions are 
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Table 4: Expected l-a upper bound on the gluino two body branching fractions at the 
example point of the AMSB model. 



Fiducial values 


c [fir 1 ] 


B(g gW°) 


B(g ^5°) 


(B(ttW),B(bbW),B(tbW)) = (0.25,0.25,0.5) 


25 


0.01 


0.02 


B(qqW) = 1 


25 


0.03 


0.02 


B(qqB) = 1 


25 


0.02 


0.06 
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enhanced. Therefore the measurement of this branching fraction provides us information 
on the scalar mass scale, which is an essential element for the Higgs mass estimation. We 
show that the measurement is reasonably possible and we can determine the stop mass 
scale. In this paper, we adopt a simple event number counting through common-used 
event cuts. Actually, more fine binning of the SUSY events can improve the discrimi- 
nation by a few factor. Or more complicated variable such as M e ^^/M e ^^ distribution 
can also help. The present method mainly relies on the event topology and not on the 
detailed of the mass spectrum. Therefore this method does not depend on the details 
of the mass spectrum. Even though there is still room for improvement, the resultant 
expected accuracy of the branching fraction is good enough to estimate the scalar mass 
scale. Although our result is based on the fast simulation, we expect that a full simulation 
can also reproduce the present result. 
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